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ABSTRACT 

Transonic Euler computations, based on unstructured grid methodology, are performed for a proposed 
High Speed Civil Transport (HSCT) configuration, designated as the Reference H configuration within the 
High Speed Research (HSR) Program. The predicted results are correlated with appropriate experimental 
wind-tunnel data for the baseline configuration with and without control surface deflections for a range of 
angle of attack at M rXj = 0.95. Good correlations between the predictions and measured data have been 
obtained for the longitudinal aerodynamic characteristics on the baseline configuration. The incremental 
effects in the longitudinal aerodynamic characteristics due to horizontal tail deflections as well as wing 
leading-edge and trailing-edge flap deflections have also been predicted reasonably well. Computational 
results and correlations with data are also presented for the lateral and directional stability characteristics 
for a range of angle of attack at a constant sideslip angle as well as a range of sideslip angles at a constant 
angle of attack. In addition, the results are presented to assess the computational method performance and 
convergence characteristics. 


SYMBOLS 

b re f reference span, inches 

c reference chord, inches 

Cd drag coefficient, Drag/(^ oc S re f) 

C Do skin friction drag at Cl = 0 

Cl lift coefficient, Lift /(go© S re f) 

Ci rolling moment coefficient, Rolling moment/ (goo S re f b re f) 

C m pitching moment coefficient, Pitching moment /(q^ S re f c) 

C n Yawing moment coefficient, Yawing moment /(^ S re f Kef) 

C p pressure coefficient, (p — Poc)/<?oo 

M a o freestream Mach number 

p static pressure 

Poo freestream static pressure 

qoo freestream dynamic pressure 

R ft Reynolds number per foot 

S re f reference area, sq. inches 

a angle of attack, degrees 

P sideslip angle, degrees 

6h t horizontal tail flap deflection angle, degrees 

5i e f inboard/outboard wing leading-edge flap deflection angle, degrees 

S te f inboard/outboard wing trailing-edge flap deflection angle, degrees 

Abbreviations 

HSCT High Speed Civil Transport 

HSR High Speed Research 

IGES Initial Graphics Exchange Specification 

LaRC Langley Research Center 
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COMPUTATIONAL OBJECTIVE AND APPROACH 

The primary objective of the computations is to provide an accuracy assessment of an unstructured grid 
Euler method known as USM3D (Ref. 1) for predicting the longitudinal and lateral/directional aerodynamic 
characteristics for the baseline configuration and the corresponding effects due to control surface deflections 
at transonic speeds. 

A computational matrix is developed to complement the configurations and the conditions tested In 
the NASA LaRC 16-Foot Transonic Tunnel. The wind tunnel test was conducted on a 1.675% model (also 
referred to as controls model) at an average Reynolds number of about 4 million per foot. Three config- 
urations are selected for the computational method assessment study. Consistent with the experimental 
test, each configuration incorporated an extended circular aftbody to represent the wind-tunnel model sting 
apparatus and appropriate boundary conditions at the nacelle’s inlet and exhaust planes to simulate the 
flow through propulsion effects. The selected configurations are: 

1) transonic cruise (i.e., Si e f = 0/10; S te f — 0/3; Sht = 0). This geometry is considered as the baseline 
transonic cruise configuration. 

2) transonic cruise with deflected horizontal tail (i.e., 5i e f = 0/10; 5 te f = 0/3; Sht = 6). 

3) supersonic cruise configuration where all the control surface deflection angles are set to zero (i.e., 
Sief = 0/0; S tef = 0/0; S ht = 0). 

The above nomenclatures are defined as: the Greek letter S represents control surface deflection angle; the 
subscripts ?e/, te/, and ht denote the leading-edge flap, trailing-edge flap, and horizontal tail, respectively. 
If applicable, the control surface deflection angles, separated by a slash, correspond to those of the inboard 
and outboard wing, respectively. For example, Si e f = 0/10 denotes that the inboard- and outboard-wing 
segmented leading-edge flap deflection angles are set to 0° and 10°, respectively. The positive control surface 
deflection angles ($i e f, S te fi and Sht ) are downward. 

The corresponding flow conditions for the selected configurations are then identified from the exper- 
imental data analyses at M c 0 = 0.95. Zero sideslip analyses are performed for all three configurations at 
a = 0°,4°, and 8° while the finite sideslip matrix composed of an a-sweep (0 o ,4 o ,8°, 11°) at /? = 3° and a 
/3-sweep (0 o ,3°,6°) at a = 4°. The selected computational matrix is depicted graphically in a chart form 
(figure 1) to illustrate the links between various configurations and the corresponding flow conditions. It 
should be noted that the transonic experimental data with finite sideslip angles were only available for the 
supersonic cruise configuration when the present computational analysis was initiated. In addition, the 
experiment was conducted to only measure the overall forces and moments. 

COMPUTATIONAL GRID 

The HSR Reference H initial surface geometry was defined in a format known as Initial Graphics 
Exchange Specification (IGES), (Ref. 2), which served as the database for all the subsequent grid generation 
processes. The GRIDTOOL (Ref. 3) program was used to read in the initial IGES geometry definition and 
the corresponding discretization into surface patches. The resulting surface patches were then read in the 
VGRID (Ref. 4) program to generate the initial surface triangulation by the advancing front method. The 
initial surface triangles were then read back into GRIDTOOL for the projection onto the initial database 
to preserve the proper surface curvature within the interior of a given patch. As a final step, the resulting 
front on the surface and the farfield boundary patches were used by VGRID to generate the volume grid 
within the computational domain along with the corresponding grid and face connectivity files. A nearfield 
and a closeup view of a typical grid, used to compute the flow for the transonic cruise configuration with 
deflected horizontal tail (i.e., 8i e f = 0/10; St e j = 0/3; Sht — 6), are shown in figures 2 and 3 from two 
different perspectives. The figures illustrate the triangular meshes on the surface, plane of symmetry, and 
the wind-tunnel sting representation. The farfield boundaries, clipped in the figures for clarity, are extended 
to about 6-longitudinal and 3-radial body lengths away from the surface. The computational grid size for 
each of the three configurations are as follows: the transonic cruise configuration had about 1.023 million 
cells and 36,000 surface triangles; the transonic cruise with deflected horizontal tail configuration had 0.822 
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million cells and 44,000 surface triangles; and the supersonic cruise configuration had 0.997 million cells and 
47,000 surface triangles. The grid for the supersonic cruise configuration was mirrored about the symmetry 
plane to generate the computational grid for all the finite sideslip runs. The resulting grid approximately 
doubled the grid size used for the zero sideslip (only one-half modeled) computations. 

ALGORITHM PERFORMANCE AND CONVERGENCE 

USM3D computations were performed on both the Numerical Aerodynamic Simulation (NAS) Cray- 
C90 located at NASA Ames, and the NASA Langley Cray-YMP. All the zero sideslip computations were 
conducted on the Cray-C90 using the implicit time integration scheme and the corresponding finite sideslip 
cases were executed on the Cray-YMP by using the explicit algorithm. Also, no flux limiters were used to 
obtain the solutions presented here. As a general strategy, solutions based on the second-order accurate flux 
difference splitting (FDS2) method were first sought; however, for certain cases only the first-order accurate 
flux difference splitting (FDSl) could be achieved. In general, the aerodynamic effects due to the difference 
in the order of solution accuracy were found to be small, particularly in terms of total forces and moments 
and will be discussed later. 

The selected computational matrix, along with the corresponding numerical method employed to con- 
verge a given solution for each configuration, are summarized in Table 1. In this table, symbols V 5 and 
‘X 1 indicate that a solution convergence ‘was’ and ‘was not’ achieved, respectively. For zero sideslip cases, 
converged solutions with FDS2 were obtained for all configurations at a = 0° and 4°, except for the super- 
sonic cruise configuration at a = 4° where only the FDSl solution could be achieved. At higher angles of 
attack (i.e., a > 4°), no solution convergence based on FDS2 could be achieved for any of the configurations 
investigated. It is speculated that the lack of solution convergence could be due to the flow separation 
occurring within the boundary layer over the round leading-edges of the inboard wing which is known to 
require Navier-Stokes simulation to resolve. As a result, the numerical solutions for a = 8° are obtained 
based on FDSl for all the configurations except for the supersonic cruise configuration for which no solu- 
tion could be achieved. The supersonic cruise configuration was the most challenging geometry to obtain 
solution convergence at the present transonic flow conditions. The solution convergence difficulties were 
identified to have been caused primarily by geometrical features particularly with respect to the transition 
region between the inboard and the outboard wing. This transition region represents a step discontinuity 
from a round to a sharp leading edge of the undeflected flap on the outboard wing. Progress has been made 
recently to implement various flux limiters in USM3D as numerical dumping mechanisms to alleviate such 
convergence difficulties. As mentioned earlier, the supersonic cruise configuration was the only geometry for 
which experimental data were available at transonic speeds with finite sideslip angle. As a result, with the 
lessons learned from the zero sideslip computations, only FDSl solutions were sought for the finite sideslip 
analysis. A typical zero sideslip case required about 500 cycles to converge and, for finite sideslip cases, 
required about 1000 cycles. The implicit time integration scheme required ~ 180 words of memory per cell 
and w 30 /x-second/cell/cycle on the Cray-C90. Similarly, the explicit scheme required ~ 46 words of 
memory per cell and ~ 6 ^-second/cell/cycle on the Cray-C90. Very little overhead CPU time, generally 
less than 5%, is found to be associated with the second-order accurate solution, compared to the first-order 
formulation. 

Typical convergence characteristics (i.e., Residuals and Lift) for solutions based on FDS2 and FDSl 
at a = 0°, j3 = 0°, for the supersonic cruise configuration are shown in Figures 4(a) and 4(b), respectively. 
A typical second-order solution is initiated automatically within USM3D after a certain level of residual 
convergence has been achieved by the first-order accurate flux formulation. For example, as shown in figure 
4(a), the FDS2 solution is initiated from the FDSl solution (after 205 iterations) where the residuals have 
already dropped 2.5 orders of magnitude. Similarly, figure 4(b) shows the residual and lift characteristics for 
the fully converged solution based on the FDSl formulation. These two particular solutions were obtained 
as part of the solution accuracy assessment between first- and second- order flux differencing schemes. A 
similar set of convergence history plots axe shown in figure 5 for typical solutions obtained for the supersonic 
cruise configuration for two different sideslip angles at a = 4°, 
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RESULTS AND DISCUSSIONS 


Surface Pressure Coefficients: 

Typical surface pressure coefficient contours computed at a = 4°, — 0.95, and (3 = 0° for the 

transonic cruise, transonic cruise with deflected horizontal tail, and the supersonic cruise configurations, 
axe shown in figures 6, 7 and 8, respectively. All surface pressure coefficients are contoured over a constant 
range from 0.6 to -1.4 with 0.067 increments (i.e., C P max = 076, C prnin = -1.4, A C p = 0.067) and shown 
from the same vantage point. The surface pressure contours computed for the transonic cruise configuration, 
shown in figure 6, serve as the baseline solution to be correlated with the results obtained for the transonic 
cruise configuration with deflected horizontal tail and the supersonic cruise configuration. The latter two 
correlations show the effects of horizontal tail deflection and the wing leading-edge and trailing-edge flap 
deflections on the computed surface pressure coefficients, respectively. Furthermore, since no experimental 
surface pressure data is available, these correlations are only intended for relative comparisons between 
computational results to show, qualitatively, the effects due to a particular control surface deflection. 

The surface pressure coefficients computed for the baseline transonic cruise configuration (figure 6) 
can be characterized by two dominant flow features, one that can be attributed to typical flow expansions 
caused by the physical geometry and the other to shock waves due to the transonic flow effects. The typical 
flow expansions on the upper surface (figure 6(a)) occur at the leading edges, particularly on the inboard 
wing, and the hinge lines associated with the deflected leading-edge and trailing-edge flaps on the outboard 
wing. The transonic flow shock waves (i.e., highlighted by a concentration of several pressure contour lines) 
on the upper surface appear in two locations. The first location is very close, nearly spanning in parallel, 
to the entire wing trailing edge, and the other location is situated approximately along the mid-chord, 
spanning the horizontal tail and extending to the fuselage and onto the vertical tail. The lower surface 
pressure contours (figure 6(b)) indicate fairly benign flow characteristics with the exception of a shock wave 
that appears in between the two nacelles and the compression rings on the wing just ahead of each nacelle 
caused by the pylons. . . _ _ 

The comparative assessments of the results shown in figures 6 and 7 indicate only very small effects 
due to the horizontal tail deflection on the computed surface pressure coefficient contours for the transonic 
cruise configuration. The most noticeable change appears to be the flip of the shock wave from the upper 
surface of the undeflected horizontal tail to the lower surface of the deflected horizontal tail. Similarly, a 
comparison of the results shown in figures 6 and 8 indicate the effects of wing leading-edge and trailing-edge 
flap deflections on the computed surface pressure coefficient contours for the transonic cruise configuration. 
As expected, the comparison indicates that the hinge line flow expansions disappear on the upper surface of 
the undeflected wing leading-edge and trailing-edge flaps. In addition, the shock waves on the wing upper 
surface near the wing trailing edge and the horizontal tail, for the supersonic cruise configuration, are not 
as crisply predicted as those of the transonic cruise configuration. This smearing of the shock waves is 
speculated to have been caused primarily by the excess dissipation generally associated with the first-order 
accurate solutions. The comparison shows very little difference on the computed lower surface pressure 
coefficient contours due to wing leading-edge and trailing-edge flap deflections. 

Forces and Moments and Correlations with Data: 

All predicted inviscid drag coefficients have been corrected with an appropriate Cd 0 to account for 
the skin friction using the experimental data obtained in the NASA LaRC 16-Foot Transonic Tunnel. This 
C&Q was determined by a linear extrapolation from the plotted C\ vs. Co curve for the supersonic cruise 
configuration as shown in figure 9. This Coo provides a minimum skin friction drag which, by definition, is 
independent of the configuration camber and twist. Through this analysis, a Coo of 0.0126 was obtained, 
which was subsequently added to all the inviscid drag predictions. 

Solution method assessment - As discussed earlier, the supersonic cruise configuration was the most 
challenging geometry for obtaining solution convergence at the present transonic Mach number. As a result, 


4 



the FDS2 solution could only be achieved for a = 0°, and the FDSl solution could only be achieved for 
a = 4°. This prompted a solution method assessment study to be performed at a = 0° for which both first- 
and second- order solutions could be obtained. The results from this study are presented in figure 10 for 
the longitudinal aerodynamic characteristics. The correlations at a = 0° indicate that the solution order of 
accuracy has only a small effect on the predicted lift and drag, with relatively more pronounced effect on 
the pitching moment estimate. 

Horizontal tail defle ction effects - The predicted and measured longitudinal aerodynamic effects due to 
the horizontal tail deflection at Mqq = 0.95 and j3 = 0° are shown in figure 11. The measured aerodynamic 
effects due to the horizontal tail deflection on both the lift and drag coefficients appear to be small. The 
pitching moment characteristics remain nearly the same over the measured range of flow conditions, and 
as expected, the horizontal tail deflection causes roughly a constant positive shift in the pitching moment 
coefficients. The overall lift and drag coefficients have been predicted reasonably well with the Euler method 
both in terms of incremental effects due to the horizontal tail deflection and magnitudes. Although the 
overall magnitudes of the pitching moment coefficients a re not predicted as well, the incremental effects due 
to the horizontal tail deflection are predicted reasonably well. 

Wing leading-edge and trailing-edge flap deflection effects - The predicted and measured longitudinal 
aerodynamic effects due to the leading-edge and trailing-edge flap deflection at M 0 0 = 0.95 and j3 = 0° 
are shown in figure 12. The measured aerodynamic effects due to the leading-edge and trailing-edge flap 
deflections appear to be small for both the lift and drag coefficients and more pronounced for the pitching 
moment characteristics particularly at ct < 10°. These aerodynamic effects due to the leading-edge and 
trailing-edge flap are reasonably predicted both in terms of overall magnitudes and trends. 

Finite sideslip stability analysis - The longitudinal aerodynamic characteristics computed for an a- 
sweep (i.e., 0 ,4°, 8°, and 11°) at 0 = 3° and = 0.95, along with the experimental measurements for 
the supersonic cruise configuration, are shown in figure 13. The correlations indicate excellent agreement 
between the computed and measured data throughout the examined range of flow conditions for both 
magnitudes and trends. 

The lateral/directional stability characteristics computed for a /?-sweep (i.e., 0°,3°, and 6°) at a = 4° 
and Mqo = 0.95, along with the experimental measurements for the supersonic cruise configuration, are 
shown in figure 14. The correlations indicate a good agreement for the pitching and rolling moments but 
not as good for the yawing moment, particularly with respect to the slope (i.e., C n0 ). 

CONCLUDING REMARKS 

The application viability of an unstructured grid Euler method (USM3D) for predicting the longitudinal 
and lateral/directional stability characteristics of the HSR Reference H configuration at a transonic Mach 
number of 0.95 have been addressed. Particular emphasis was placed on assessing the ability of the method 
to provide consistent and accurate solutions of forces and moments for the baseline configuration with 
different control surface deflection settings for the horizontal tail as well as for the leading-edge and trailing- 
edge flaps. The analysis revealed that a consistent set of solutions based on the first-order accuracy was 
achievable for a wider range of applicability (in terms of configurations and flow conditions) where the 
second-order accurate solutions could not be obtained. The solution method assessment analysis indicated 
that the first-order accurate solutions provide comparable results to those of the second-order, particularly 
for the aggregate forces and moment characteristics. The computed longitudinal aerodynamic effects due to 
control surface deflections and sideslip angles agree reasonably well with the measured data both in terms 
of magnitude and trends. Finally, the lateral/directional stability analysis revealed fairly good agreement 
between the predictions and measured data for the pitching and rolling moment characteristics but not as 
good for the yawing moment. It is believed that the present unstructured grid Euler method is certainly a 
viable engineering tool that can be utilized with confidence in the early aerodynamic design and analysis of 
the proposed high speed civil transport configuration. 
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Table 1.- Computational matrix - method of convergence. 



Figure 1.- Selected computational matrix - geometries and conditions. 
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(b) Closeup. 


Figure 2.- Typical computational grid - Upper surface view. S te f = 0/10, 6 te j = 0/3, 8 ht - 6. 
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(b) Closeup. 


Figure 3.- Typical computational grid - Lower surface view. Si e j = 0/10, $tef — 0/3, Sht — 6. 
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Figure 5.- Typical FDSl convergence characteristics - 8 lef = 0/0, 8 te f = 0/0, 6ht = 0; Af™ = 0.95, a = 0 
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Figure 9.* Cdo determination strategy - Supersonic cruise geometry, — 0.95 
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Figure 10.- Solution method accuracy assessment - Supersonic cruise geometry, M a Q — 0.95,/? = 0°. 
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Figure 11.- Horizontal tail deflection effects and predictions - Transonic cruise geometry, Moo = 0.95, (3 = 0°. 



Figure 12.- Leading/trailing-edge flap deflection effects and predictions - Transonic cruise geometry, — 0.95, (3 = 
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Figure 13.- Longitudinal aerodynamic characteristics and predictions - Supersonic cruise geometry, M co = 0.95,/? = 3°. 



Meg Meg 

igure 14.- Lateral/directional stability characteristics and predictions - Supersonic cruise geometry, M <*> = 0.95, a = 4°. 
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